Adaptive Control of Scalar Plants in the Presence of 

Unmodeled Dynamics 

Heather S. Hussain, Megumi M. Matsutani, Anuradha M. Annaswamy and Eugene Lavretsky 



o 



m 

00 

o 



> 

o 
o 

(N 
O 
cn 



X 
S3 



Abstract — Robust adaptive control of scalar plants in the pres- 
ence of unmodeled dynamics is established in this paper. It is 
shown that implementation of a projection algorithm with standard 
adaptive control of a scalar plant ensures global boundedness of 
the overall adaptive system for a class of unmodeled dynamics. 



I. Introduction 

Following the Rohrs counterexample in |i9], several robust 
adaptive control solutions were suggested in the '80s and '90s 
(see, for example, |7| and [2|), including specific responses to 
the counterexample (see for example [1], [3J, [8J, [2], and [6J). 
Most of these were qualitative, or local, and often involved 
properties of persistent excitation of the reference input. In 
this paper, we show that for a class of unmodeled dynamics 
including the one in jO), adaptive control of a scalar plant with 
global boundedness can be established for any reference input. 
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Fig. 1. Adaptive control in the presence of unmodeled dynamics 



II. The Problem Statement: Scalar Plant 

The problem we address in this paper is the adaptive control 
of a first-order plant 

Xp{t) = apXpit) + v{t) (1) 

where ap is an unknown parameter It is assumed that \ap\ < a, 
where a is a known positive constant. The unmodeled dynamics 
are unknown and defined as 



v{t) = Cr,'^Xr,{t) 



(2) 
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Xjj{t) is the state vector, and u{t) is the control input. The goal 
is to design the control input such that Xp{t) follows Xm{t) 
which is specified by the reference model 



X„i{t) = araXm{t) + r{t) 



(4) 



where a,„ < 0, and r{t) is the reference input. The adaptive 
controller we propose is a standard adaptive control input given 
by (see figure [T]i 



u(t) = e{t)xp{t)+r{t) 



(5) 



where the parameter 9{t) is updated using a projection algorithm 
given by 



where 



^(t) =7Proj(0(t),-Xp(t)e(t)), 7>0 

e(t) = Xp{t) - Xm{t) 



Proj(6l,y) 
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y [6 e Qa, y0 > 0] 



otherwise 



(6) 



(7) 



(8) 



no = {0eR'\ -Ca. < < CaJ 
= {0 e I —9„iax < < Omax} 
— ^l\^0 

with positive constants 0'^^^ and 9„iax given by 

Knax >a+\am\ 
Qmax — Omax 



(9) 

(10) 
(11) 



Lemma 1. Consider the Adaptive Law in dSJl with Projection 
Algorithm in dS]) to Mlh Then, 



ma)\\ < 



\m\\ < 



yt > ta. 



(12) 



Hence, the projection algorithm guarantees the boundedness of 
the parameter 6{t) independent of the system dynamics. We refer 
the reader to [t4J for the proof of Lemma [7] 



III. Choice of Projection Parameters 
The projection algorithm in (|8) is specified by two parameters 
^'max ™d 9max- Equation ( fTOl i provides the condition for O'^^x- 
To determine Eq in ( fTTT i. the following discussions are needed: 
We consider the linear time-invariant system specified by ([T]i, 
dU, and (|5]i, with the parameter 9{t) fixed as 

0it) = -0max, yt > ta. (13) 

The closed loop transfer function from r{t) to Xp(t) is given by 



Gc{s) 



9c(s) 



where Gc(s) is defined using G,j{s) in (O as 

9c(s) = <7,,(s)(s - Op) + 6'maa;Pi,(s). 

From ( [Tol l and (fTTT i. it follows that 



^7nax 



< 0, V|ap| < 



(14) 

(15) 
(16) 

(17) 



Therefore it follows that there exists a class of unmodeled 
dynamics (c^, A^j, fe,,) such that qds) has roots in C^, the left- 
half of the complex plane. It is this class that is of interest in 
this paper. 

Definition 1. The triple (c^,A^,6^) is said to belong to 
Sri{a,9max) if Pri{s) and qrf{s) in Iil5\l are such that the roots 
of qc{s) in ( 1761 ) lie in for all \ap\ < a. 

Let's demonstrate 5^ (a. 9 max) with an example. Consider the 
class of unmodeled dynamics of the form 

Gni^) = 2^of^ ^ 2 (18) 

where C > and a;„ > 0. From ([Till, GSll, and ([T8]l, the closed- 
loop dynamics from r to Xp is given by 

GM - ^ (19) 



where 



9c(s) 

qc{s) = + ais^ + a2S + 03 
ai = (2Cw„ - Qp) 
a2 = {ll>1 ~ 2apCuj„) 

0,3 



(20) 



For the roots of qds) in ( |20l ) to lie in C^, the following 
conditions are neccessary and sufficient for all \ap\ < a: 

(A-i) ttp < min(2Ca;„, —y ] 

(A-ii) 9rnax > O-p 



(A-iii) 
If 

where 



< 



2C' 



~r = (2Cw„-ap)( 



^ _ 2(ap , 



(21) 
(22) 



then conditions |(A-ii)| and |(A-iii)| hold. 

Hence, any class of unmodeled dynamics (c^,A^,6^) in 
(fTsT l satisfying condition |(A-i)| belongs to Sjj{a,6max)- It can 
be easily shown that the unmodeled dynamics and the plant 
discussed in the infamous Rohrs counterexample |l9] satisfies 
conditions |(A-i)| to |(A-iii)| above for some dmax- 

We now discuss the choice of Eq. Consider the class of 
unmodeled dynamics Sjj{a,6max) in Definition [1] Since the 
closed loop system specified by ([T]i, dU, (|5]l, and (fT3] l is stable, 
it follows that there exists a Lyapunov function 



with a time derivative 



V = -x^Qx 



(23) 



(24) 



where ^ ] . P is the solution to the Lyapunov equation 



A^P + PA = -Q <0 



where 



A^ 



Ar 



(25) 
(26) 

satisfies d^TT l. 



since A is Hurwitz. The latter is true since 0„ 
We define two sets Vlu C VLa and VLi C VIa as 

^u = {9eR^\ -Bmax + eo < < -9'max} (27) 

= e I —9raax < 6* < —9max + Co} (28) 

where 

eo=c£o, ce(0,l). (29) 

We now consider the linear time-varying system specified by 
©, (Ell, and dSll, with 9{t) e r2„ U rj;. It follows from (EJ and 
(fT2] | that 



e{t) = -9„,ax + e{t), \f9{t)enuuni 
d{t) = -9,nax + m, "^9(1) eni 



where 



E{t) e[o,£o), m&[0,^o]- 

Therefore, the closed-loop system is given by 

X ^ Ax + A^{t)x + br, \/9{t)eni 

where 

Mt) 

If we choose V — - 

V < -Ao^„ 



0' 




■q" 
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(30) 
(31) 

(32) 

(33) 
(34) 



—x^Qx with P as in (l25T l, we obtain 

2 + 2Ap_fceol|sf + 2Ap_ 



where 



That is. 



Aq„„„ ^min|K(A,(Q))| 
Ap„_ ^max|3?(A,(P))| 
|6»7ll < k, r,nax = max|r(t)|. 

F<0 if||x||>a;o 



'^max 1 1 *^ 1 1 

(35) 



(36) 
(37) 

(38) 
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where 



2A, 



Xo 



A = A 



2\f 



Mo- 



rn 

(40) 



In summary, the closed-loop system has bounded solutions for 
all d{t) G ill with ||a;(t)|| < a;o if (c^,A^,6,,) is such that 
(B-i) qc{s) has roots in for all \ap\ < a, and 
(B-ii) ^0 < eo , where 

We introduce the following definition: 

Definition 2. The triple (cj,, is said to belong to 

Sri{a,9max-,£,o) if Conditions \( B-i)\ \( B-ii)\ and \(B-iiTj\ above are 
satisfied. 

IV. Main Result 

Theorem 1. Let z{t) = [e{t) 0{t)]'^. The closed-loop adaptive 
system given by (Ii}-(1Z3 has globally bounded solutions for all 

e{ta) G if {Cr,,Ar„br,) G Sr,{a, Oynax , 

Definition 3. We define the region A and the boundary regions 
B and B as follows 

B = {zeR^\e:^^^<9<0^a4 

{z eW' \e eflo} (41) 

B = {z G I -9max <0< -0'^^^} 

Definition 4. We divide the boundary region B into two regions 
as follows: 

Bu ^{z eM.^ \ 9 e n.u} 
Bl = {zgR2 \ 0eni} 

with B ^ Bu D Bl- 

z(t) = [e{t) e(t)f 



(42) 
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Fig. 2. Definition of regions in I4H and j42t . and phase.s I-III 
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By combining the adaptive law in Q and and boundary 
region definitions in dTTl l. we obtain 



6*2 



-^exp if z G (S UB), -eXpO > 



otherwise 



(45) 

Hence, the projection algorithm guarantees the boundedness of 
the parameter d{t) independent of the system dynamics|4|. It 
follows that Theorem [1] is proved if the global boundedness of 
e{t) is demonstrated. This is achieved in four phases by studying 
the trajectory of z{t) for all t > ta- This methodology was 
originally proposed in ||5| for adaptive control in the presence 
of time delay. 

We begin with suitably chosen finite constants e and 6 such 
that e — 6 > 0. The trajectory then has only two possibilities 
either (i) \e{t)\ < e — S for all t > ta, or (ii) there exists a time 
ta at which |e(ta)| = e — S. The global boundedness of e{t) 
is immediate in case (i). We therefore assume there exists a ta 
where case (ii) holds. 

(I) Entering the Boundary Region: We start with |e(ta)| — 
e — 5. We then show that the trajectory enters the boundary 
region _B at G {ta, ta + ^Tb), and B^ at tc > t), where 
ATb and t^ are finite. 
(II) In the Boundary Region, B^: When the trajectory enters 
B, the parameter is in the boundary of the projection 
algorithm; e is shown to be bounded in Bl by making 
use of the stability property of the underlying linear time- 
varying system. For t > tc, the trajectory has only two 
possibilities: either (i) z stays in Bl for all t > tc, or (ii) 
z reenters Bjj at some td > tc where |e(fd)| < Xm- 

(III) In the Boundary Region, Bjj'- For t > td, the trajectory 
has three possibilities: either (i) z reenters A at t — tc, 
(ii) z stays in Bjj for all t > td, or (iii) z reenters Bl at 
tf G {td,td + ATbl) where ATg^ is finite. 

(IV) Return to Phase I or Phase IT. If case (i) from Phase III 
holds, then the trajectory has only two possibilities: either 
\e{t)\ < e — S, \/t > tc which proves Theorem [T] or there 
exists atg > tc such that |e(fg)| = e — 5 in which case the 
conditions of Proposition [T| are satisfied with ta replaced 
by tg, and Phases I through III are repeated for t > tg. If 
case (ii) from Phase III holds, then the boundedness of e 
is established for all t > td- If case (iii) holds, then Phases 
II and III are repeated for t > tf. In all cases, e remains 
bounded throughout. 



Proof of Theorem [7] The closed-loop adaptive system has error 
dynamics in (|7| equivalent tcQ 



e = QmB + Sxp + r] 



(43) 



where 



u = u — a , a = Gm ^ CLp, rj = V — u. 

'For ease of exposition, we suppress tlie argument "t" in wliat follows. 



A. Phase I: Entering the Boundary Region 

We start with |e(ta)| — e — 6. From (|44] |. it is easy to see 



that 



\t]\ < (fc,, + 1)6' 

7nax \ \^\ < "^ni 

(44) where A:^ = ||G^(s)|| and 



max \xm{t)\. 
t>ta 



(46) 



(47) 
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We define e as which implies that 

e = max{eo, ei} (48) 



\e{t)\<e, yte[ta,ta + AT] (64) 
and the proof of Proposition [lUi)| is complete. □ 



where 

eo = \xp{t + 2S (49) 

e^ = l (cbQ + Jc^bl + 4c6i) (50) Proof of PropositionUH^ We note from (ESll that 



2 



■fuu AKA^ A- ^ HR7T,A^^n- ^ ^rn 1 |e(i)| > |e(t„)| - + 6i)AT, G [t,,, ia + AT]. (65) 
with Oo and Oi denned m (153JI and (154)1, d e (0, Xm), a G (0, eoj, 

c in dllll, and Sin^-g |e(i„)| = e - 5, (|65l) can be simplified as 



^. (51) 



|e(i)| > e-2(5, e [ta, ia + AT]. (66) 
Since e > |xp(tQ)| + a;™ + 2(5 and (5 < Xm, it follows that 



(57 

Phase I is completed by proving the following Proposition: 
Proposition 1. Let z{ta) G A with \e{ta)\ = e — 6 where e is e — 26 > Xm- (67) 



This in turn implies that 9{t) is negative for all i e [ta , ta + AT] 
with 



given in ( 14^1 ) and 5 G {0,Xm)- Then 

(i) \e{t)\ <e, Vte [ta,ta + AT] 

(ii) z{tc) G Bl for some U G (ta,ta + AT) 
where 9{ta) - e{ta + At) > -f{e - 2S){e - 26 - x,n)At (68) 

AT = (52) 

boe + bl for all At G [0, AT]. Defining, 

foo = |am| + (fc„ + 2)(?„.,, + |r| (53) _ ^^-""^ -eo + a ^^^^ 

fei = ((fc„ + 2)d„,a. + \e*\)x„, + (fc,, + 2) w (54) 7(e - 25)(e - 2,5 - S™) 

Proof of Proposition Um From (|43ll and it follows that it follows that z{tb) enters B at 4 e {ta,ta + ATg) if ATs < 

AT. 

|e(t)| < 6oe' + bl, Vi G [^a,ta + AT] (55) We now show that z{tc) enters at t^ < t^ + AT'g for 

^jjgj.g some AT^ > ATg. It can first be proven that 

= ^^^^ |Proj(0,y)|>c|z/|, yzeBu. (70) 

We will show below that e' — e which proves Proposition [TJi). xhen from (|45]l 
We have that for all At G [0, AT], 

\e{ta + At)\<\e{t,)\+ max \e{t)\^T (57) - m > ic{e - 25){e - 25 - x^), yt e Tb,. (71) 

te[ta,ta+AT] 

: where Tb„ is defined as 

< {e ^ S) + ib„e' + bl) (58) 

+ Tbu ■■ {t I z{t) G S[/ and i G [ta,ta + AT]}. (72) 
from (|55j, ( |56] l. the definition of AT, and the choice of \e{ta)\- 

From (|56ll, with some algebraic manipulations, ^ can be Since the distance the trajectory can travel in Bl is bounded by 

rewritten as Co^ the maximum time z{t) spends in Bjj can be derived from 

e' < e + (g/ _ g) (59) GB, and we obtain 

ooe + 6i 

which can be simplified as ATb^ = — jz , , — j ; — r. (73) 



(e' - e) (1 - &oAT) < 0. (60) 



7c(e — 2(5) (e — 2(5 — Xm) 
This implies that z{tc) enters Bl at tc G {ta,ta + ATg) where 



Since 6 < Xm, from the definition of bo and bi, it can be shown 
that AT^j - ATb + ATs^ (74) 

(1 - 6oAT) > 0. (61) „ 

if AT^ < AT, since then dlB is satisfied for all t G (tb,tc]- 
Therefore from dgO), it follows that From the choice of e in (|48]l, we have that 



e' - e < 0. (62) 1 , 

e > - ( c6o + .Jc^bl + Acbi ) . (75) 
From the definition of e' in (|56] |. it follows that only the equality ^ 



in (I62ll can hold. Hence 



Using algebraic manipulations, it can be shown that (175b implies 
e{ta + At)\< e, VAi G [0, AT] (63) that AT^ < AT. This proves Proposition [In)] □ 
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B. Phase II: In the Boundary Region, B 

When the trajectory enters Bl, the parameter is in the 
boundary of the projection algorithm with thickness i^o; e(t) 
is shown to be bounded by making use of the underlying linear 
time-varying system in (|33] | and (|34] |. 

Let z{t) e Bl for t e [tc,td). That is, 0{t) = -0max + 
^{t) for t e [tc,td) with {{t) satisfying ^ and Since 
{cr,,A,^,br^) e Sr^{a,9max,£.o), from (l38T l. it follows that 



Mt)\\<xo, ^teTB^ 

where Tb^ is defined as 

Tsr, ■■ {t I z{t)eBL}. 



(76) 



(77) 



Since |e(i)| < \xp{t)\ + Xm for all t G Tbj^ and x — [xp x'^]'^, 
this implies 

|e(i)|<e2, yt£{tc.,td) (78) 



where 



62 = 2^0 + a;„ 



(79) 



which proves boundedness of e in S^. 

We have so far shown that if the trajectory begins in A at 
t ~ ta, it will enter the region Bj^ att — t^, where < ia+AT, 
and AT is finite. For t > tc, there are only two possibilities 
either (i) z stays in Bl for all t > tc, or (ii) z reenters Bjj at 
t — td for some td > tc- If (i) holds, it implies that ( iTST i holds 
with td — oo, proving Theorem [T] The following Proposition 
addresses case (ii): 

Proposition 2. Let z{t) e Bl for t e [tc,td) and z{td) G Bu 

for some td > tc- Then 



\e{td)\ < Xr, 



(80) 



Proof. Since z{t) e Bl for t e [ic7^d) and ^(td) G ^S^/ for 
some td > tc, from ( |42] |. it follows that for any Ai^j G (0, — 

tc], 

9{td — Atd) < —0,nax + £.0, &{td) > —Qmax + £.0- (81) 

This implies that 0{td) > which in turn implies 

\e{td)\<x„,. (82) 
which proves Proposition |2] □ 
We note from A79\l that 

|e(i)|<e2, yte{tc,td] (83) 
which proves boundedness of e in Phase II. 

C. Phase III: In the Boundary Region, Bjj 

The boundedness of e has been established thus far for all 
t G [^ajirf]- For t > td, there are three cases to consider: either 
(i) z reenters A ai t — tc for some tc > td, (ii) z remains in 
Bu for all t > td, or (iii) z reenters Bl at t/ G {td,td + ATb^ ) 
with ATb^ given by (l73T l. 

We address case (i) in the following Proposition. 



Proposition 3. Let z{t) G Bu for t G [td,tc) and z{tc) G A 
for some tc > td. Then 



|e(t)|<a;„, VtG(trf,ie] 



(84) 



Proof. Since z{t) G Bu for t G [td, tc) and z{tc) G A for some 
tc > td, from dTO . it follows that for any Ate G (0,ie - id], 

- At,) < 0(t,)>-Ca.- (85) 

This implies that d{t) is positive, and we obtain 

\e{t)\<xm, yte{td,tc] (86) 

which proves Proposition |3] □ 

We now address case (ii) and (iii). 

We consider suitably chosen finite constants 63 and S such 
that 63 — (5 > 0, and 



where 



63 = max{e2, 63} 



62 = 2Xm + 2(5 

63 = ^ f C260 + \/c2&0 +4C2&1 



and 



C2 



(1 - c)£o 

(57C 



From 



and the definition of 63, it follows that 

|e(td)| < e^-S. 



(87) 

(88) 
(89) 

(90) 
(91) 



If e{t) grows without bound, it implies that there exists > tc 
such that 

\e{t'c,)\ ^ es ~ S. (92) 



Hence, 



|e(t)| < 63 -(5, yte[td,t'a). 



(93) 



We show below that if such a t^ exists, then z(t) must enter 
Bl at t — tf, for some finite tf > t'^. 

Proposition 4. Let z{t) G -Bjy /or all t G [td,tf), and 3t'^ G 
(id, i/) SMc/z that \e{t'j)\ —63 — 6 where 63 is given in l[87\l and 
S G (0, Xm)- Then 

(i) \eit)\<e3, Vte[t'a,t'^ + AT'] 

(ii) z{tf) G Bl for some tf G [t'^, t'^ + AT') 
where 

AT' = , _ . (94) 

6063 + bl 

Proof. We note that Proposition (|4|i is identical to Proposition 
[T]with ta replaced by t'^, e replaced with 63, and z{t'^) G Bu 
which implies ATb = 0. Using an identical procedure, we can 
prove both Proposition |4 |ti)| and Proposition |4 Uii)| □ 

We note that if case (ii) holds, it implies that ( |93] | holds for 
t'^ = 00, which implies that e{t) is globally bounded. 

In summary, in Phase III, we conclude that if z enters Bu at 

t = td, 

(i) z enters A at t = tc with \e(t)\ < Xm for all t G [td,tc]. 
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(ii) z remains in Bu for t > td with \e{t)\ < 63—8 for all 
t > td, or 

(iii) z enters at t — tf for tj > with \e{t)\ < 63 for all 

t e [td,tf]. 

Therefore, either Phases I and II, or Phases I, II, and III, 
can be repeated endlessly but with \e{t)\ remaining bounded 
throughout. This is stated in the next section. 

D. Phase IV: Return to Phase I or Phase II 

If Proposition [3] is satisfied, then the trajectory has exited the 
boundary region and entered Region A. Therefore, \e{t)\ < e — S 
for all t > te, in which case the boundedness of e is established, 
proving Theorem[T] or there exists a tg > such that |e(ig)| = 
e — 6. The latter implies that the conditions of Proposition [T] are 
satisfied with ta replaced by tg. Therefore, Phases I through III 
are repeated for t > tg. 

If Proposition!?] is satisfied instead, then z has reentered Bl, 
in which case Phases II and III are repeated for t > tf. 

By combining (gH) from Phase I, (|78]l, and ^ from 
Phase II, and ( f86b and ( f87b from Phase III, we obtain 

|e(i)| < max{e, 62,63}, 'it>ta 

proving Theorem [T] 

V. Numerical Example 

In this section we demonstrate using the counterexample in 
||9l as to how the main result in this paper can be used to obtain 
robust adaptive control in the presence of unmodeled dynamics. 
We consider the nominal first order stable plant0 

x^) = j^^[u(t)] (96) 

in the presence of highly damped second order unmodeled 
dynamics, described by (fTsT l with 



(95) 
□ 



C = 0.9912, ijjn = 15.1327 
and a reference model 



Xm (^) 



(97) 



(98) 



(s + 3)' 

The adaptive controller is chosen as in (ISt-lfTTTi with suitably 
chosen 6„iax and Eq — 0.1f?max- The control problem differs 
slightly from that shown in Fig. [T]and requires gain compensa- 
tion on the reference input. 

That is, the plant and reference model differ from dl) and (3) 
such that 

Xp{t) — apXp{t) + kpv{t) (99) 
(t) + k^r{t) (100) 

where 

flp = -1, kp = 2, am = -3, km = 3. (101) 
The control input is then 

u{t) = e{t)xp{t) + krr{t) (102) 

in what follows is a differential operator d/dt and not the Laplace variable. 



where kr = '=™/fcp — 1.5 so as to match the closed-loop adaptive 
system when no unmodeled dynamics are present (G,j{s) = 1). 

We now show that ( fTsT i with ( |97| l corresponds to 
Sjj{ap, 9max, £q) for suitably chosen 9,nax and Eq. When d{t) — 
—Smax in (IIO2I 1. the closed-loop adaptive system given by ( |99] l, 
(llOOI l, ( llOll l. and ( I102l i has a transfer function from r to Xp of 
the form 

GM = , , ^J^^^^^,^ , . (103) 



5 + 31s2 + 259s + 229 

In addition to condition |(A-i)[ the following conditions are 
neccessary and sufficient for the poles of Gc(s) in ( |103l l to 
lie in C^, which are slightly modified versions of |(A-ii)| and 
|(A-iii)| due to the presence of kp and km'- 

(A-ii)b Qm, 

(A-iii)b Bm, 
Therefore, if 0„ 



> 



< 



2Cw, 



is such that 



< 



< -r(ap 

hip 



(104) 



with Q* = 35.06 then conditions [(A-ii)b| and |(A-iii)b| hold. Since 
ap — —1 and kp = 2, if we choose dmax = 16.7, (1104b is 
satisfied. With C and w„ in (|97| |. |(A-i)| is satisfied as well. 

We now demonstrate the choice of eo- Since the closed-loop 
system in ( 11031 ) is stable for satisfying |(B-iii)[ a Lyapunov 
function is chosen as in ( l23T l. It follows from dZST l and ( |36] | that 
Q and P are such that Ag^^^^^ = 1 and Xp^^^ — 47773.6. Since 
imil = 229 from (fTSl l and ( |97] |. we choose ^0 using |(B-iii)| such 
that 

^0=4.57-10"^ (105) 



Condition |(B-ii)| implies that any Eq such that < Eq < Omax 
suffices, with the actual value determined between the trade-off 
between adaptation and numerical accuracy. In the numerical 
simulations we report below, we chose eo = O-Winax- 

In summary, 9max — 16.7, ^0 as in ( 11051 ), and eq = 1.7, 
ensures that the triple (c^, A,,, fe,,) belongs to Srj{ap, 9 max, ^o)- 
With these choices, the adaptive controller in (l45t and (11021 ) 
guarantees globally bounded solutions for any initial conditions 
Xp(Q) and 6'(0) with ||6l(0)|| < 6*,, 
dynamics in ( fTSl l and WT\ - 



for the Rohrs unmodeled 



A. Simulation Studies 

In this section, we carry out numerical studies of the adaptive 
system defined by the plant in ( |99] l in the presence of unmodeled 
dynamics in ( fTSl l and ( |97] i with the reference model in (1 100b . the 
controller in (I102t . and the adaptive law in (l45t with 9max = 
16.7 and £0 = 1-7. The resulting plant output, Xp, reference 
model output, Xm, error, e, and 9 are illustrated in Fig. |3]for 
the reference input 



r{t) =0.3 + 1.85sin(16.1t) 



(106) 



and initial conditions Xp{Q) — and 9{{)) — —0.65. It was 
observed that all of these quantities became unstable when the 
projection bound in (fTTb was removed. It is interesting to note 
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Fig. 3. Simulation of stable plant in the presence of highly damped unmodeled 
dynamics with adaptive law in ^45). 



that in this case, only Phases I and II discussed in Section [TVl 
occurred, with Phase I lasting from t = to t = 1377.5s and 
Phase II for all t > 1377.5s. This clearly validates the main 
result of this paper reported in Theorem 1. In what follows, we 
carry out a more detailed study of this adaptive system, by only 
changing the reference input. As the numerical simulations will 
show, the behavior of the adaptive system, in terms of which 
of the four phases reported in Section |IV] occur, is directly 
dependent on the nature of the reference input. Four different 
choices of the reference input are made, and the corresponding 
behavior are described. 

(i) r{t) = 0.3 + 2.0 sin(8t): The error, e, and parameter, 6, 
corresponding to this reference input are shown in Fig. 5] 
We observe immediately that |e(t)| < 1 for all i > 0. As a 
result, the trajectory never enters B, eliminating the need 
for Phases II, III, or IV. Hence, no projection is required 
in this case. 
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Fig. 4. Simulation results with reference input in|(i)| 

(ii) r{t): A pulse for the first one second. That is. 



rit) 



12 0<t<ls 
t>ls 



(107) 



The corresponding trajectories are shown in Fig. which 
illustrate that Phase I occurs for < t < 0.9s, and Phase 



II for 0.9s < t < 1.0s. The trajectory exits the boundary 
region at = 1.0s, demonstrating Phase III. Phase I 
is repeated, and the trajectory reenters S at = 1.3s, 
demonstrating Phase IV. The trajectory then settles in B 
for all t > 1.3s. 
II/ 




Fig. 5. Simulation results with reference input in |(ii)| 

(iii) r{t) — 10 Vt: Fig. |6] illustrates the corresponding limit cy- 
cle behavior of the trajectory. We observe that the trajectory 
first enters _B at if, = 1.80s. Phase II then occurs for 
1.80s < t < 9.82s. Phase III occurs for te = 9.82s, and 
then Phase I is repeated with the trajectory reentering B 
at tb = 9.84s. Phases I through III are repeated for all 
t > 9.84, demonstrating Phase IV, a limit cycle behavior 
The points at which the trajectory enters B (i.e. Phase II) 
are shown in orange, and the points at which the trajectory 
exits B (i.e. Phase III) are shown in purple. 
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Fig. 6. Simulation results with reference input in |(iii)| 

(iv) r{t) ~ 2 sin(ti;oi), with wq undergoing a continuous sweep 
from 16.1 rad/s to 2 rad/s over a sixty second interval. The 
resulting trajectories are shown in phase-plane form in Fig. 
ItJ), for six different initial conditions labeled 1 through 6. It 
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is observed that the trajectory behaves differently for each 
initial condition. Initial conditions 1 through 3 resulted 
in trajectories that remained in Phase II for all t > tb. 
Initial condition 4 led to a trajectory with a finite number 
of occurrences of Phases I through III and finally settled 
in Region A. Initial conditions 5 and 6 stayed in Region 
A for all t > 0. All steady state values are labeled as 1/ 
through 6/. 

lOf 
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Fig. 7. Simulation results with reference input in |(iv)| 



VI. Summary 

In this paper, robust adaptive control of scalar plants in the 
presence of unmodeled dynamics is investigated. It is shown 
through analytic methods and simulation results that implemen- 
tation of a projection algorithm in standard adaptive control law 
achieves global boundedness of the overall adaptive system for 
a class of unmodeled dynamics. The restrictions on the class of 
unmodeled dynamics and the projection bounds are explicitly 
calculated and demonstrated using the Rohrs counterexample. 
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